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ABSTRACT 

In this research paper, the Natural Decomposition Method (NDM) is implemented for solving the linear and 
nonlinear Klein Gordon equations. 

The method which is based on the Natural Transform Method (NTM) and the Adomian Decomposition Method 
(ADM) is use to obtain exact solutions of three modelling problems from Mathematical Physics. 

The results obtained are in agreement with existing solutions obtained by other methods and demonstrate the 
simplicity and efficiency of the NDM. 

KEYWORDS: Adomian Decomposition Method, Klein Gordon Equations, Natural Transform, Sumudu Transform, 
Laplace Transform 

1. INTRODUCTION 

The Klein Gordon Equation is considered one of the most important Mathematical models in quantum field 
theory, with appearances in relativistic physics, nonlinear optics and plasma physics. It arises in physics in linear and 
nonlinear forms and it is useful in describing disperse wave phenomenon [1]. 

We consider the Klein Gordon Equation 

u tt (x,t)- u vc (x,t) + au(x,t)+ F(u(x,t))= h(x,t) (1.1) 

subject to the initial conditions 

u{x,0) = f{x) ; u,{x,0)=g{x) 

Where U is a function of X and t, a is a constant, h(x,t) is a known analytic function and F{u (*. 0 ) is a 
nonlinear function of u{x,t ) 

There are many integral transform methods [3-7] existing in the literature to solve PDEs, ODEs and integral 
equations. Many numerical methods were developed recently for solving Klein Gordon Equations such as Reduced 
Differential Transform Method (RDTM) [2], Adomian Decomposition Method (ADM) [8,9]and Variational Iteration 
Method (VIM) [10,11]. 

In this paper, the following Klein Gordon Equations were solved: 

First, consider the homogenous Klein Gordon Equation [10] 
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U n ~ U ,,~ U =^ 

Subject to the initial conditions 
u(x,0) = 1 + sin a ; u,(x, 0) = 0 

Where U = U\X,t) is a function of the variables X and t . 

Secondly, the in homogenous nonlinear Klein Gordon Equation [2] 
U„ — U xx +U 2 =-ACOs(f)+ X 2 COS 2 (t) 
subject to the initial conditions 

u(x,0) = x \ u,(x, 0 ) = 0 

and lastly, the nonlinear non-homogenous Klein Gordon equation [1] 
u tt~ u xx + u — - 2 1 +xt 


(1.3) 


(1.4) 


(1.5) 


( 1 . 6 ) 


(1.7) 


Subject to the initial conditions 

u{x,0) = u t (x,0) = 0 (18) 

The structure of the paper is organized as follows: In section 2, basic idea of the Natural Transform Method is 
discussed, section 3 give definitions and properties of the N-Transform. In section 4, the methodology of the NDM is 
explained and in section 5, the NDM is applied to solve three test examples in order to show its simplicity and efficiency. 
Section 6 is the conclusion. 


2. BASIC IDEA OF THE NATURAL TRANSFORM METHOD 


Here we discuss some preliminaries about the nature of the Natural Transform Method (NTM). 

Consider a function f{t\ t £ (— oo ? oo), then the general Integral transform is defined as follows [7, 12]: 


3[/M](s)= \k{s,t)f{t)dt 


( 2 . 1 ) 


where k(s,t) represent the kernel of the transform, S is the real (complex) number which is independent of t . 

Note that when k(s,t )is e '' , tJ (st) and t ' 1 {st) , then Equation (2.1) gives, respectively, Laplace Transform, 
Hankel Transform and Mellin Transform. 

Now, for f(t), t G ( — oo ? oo) consider the Integral transforms defined by 
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3[/(? )](«)= [k(t)f(ut)dt 

( 2 . 2 ) 


and 3[/ (?)](.S',m)= j k(s,t)f(ut)dt 


(2.3) 


Note that: 


• when k(t) = e ' , Equation (2.2) gives the Integral Sumudu transform, where parameter S is replaced by It . 
Moreover, for any value of n , the generalized Laplace and Sumudu transform are respectively defined by [7,12]: 


Af{t)\=F(s) = S n ]e- s " +1 'f( S n t)dt 

0 

and S[f{t)] = G(u) = «"j e- u "'f(tu" +l )dt 

o 

• when Tl = 0, Equation (2.4) and Equation (2.5) are the Laplace and Sumudu transform respectively. 

3. DEFINITIONS AND PROPERTIES OF THE N-TRANSFORM 


(2.4) 


(2.5) 


The natural transform of the function f(t\ t G (— 00 , 00 ) is defined by [7, 12]: 
N[/(f)] = R{s, u)= J e~ st f(ut)dt ; s,ue (-00,00) 


(3.1) 


where N[/ (f)] is the natural transformation of the time function /(f) and the variable S and it are the natural 
transform variables. 

Note: 

• Equation (3.1) can be written in the form [7,12]: 


N[/( ? )]= l e s 'f{ut)dt ; s,ue (-00,00) 


U '-KJ 

je~ s, f(ut)dt ; s,u e (-°°, 0 ) + je~ s 'f(ut)dt ;s,ue( 0 ,°°) 


J Lo 

= N~[/(f)]+N + [/(f)] = N[/(f)H(-f)]+ N[/(f)H(f)] = R~ (s ,u) + R + (s ,u) 
where h() is the Heaviside function. 
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If the function f (f) H (f ) is defined on the positive real axis, with t £ (0,°°) and in the set 

M 

A = \ fit) 3 M ,T l ,T 2 > 0, s.t\f(t )| < Me Ti , where te (-l) ; x[0,°°), j = 1,2 


> , then we define the 


Natural transform (N-Transform) as [3,4]: 


N[/(r)H(r)] = N + [/(f)] = 7? + (5,i<) = °je s, f(ut)dt ; s,ue{ 0,°o) (3.2) 

o 

• If U = 1, Equation (3.2) can be reduced to the Sumudu transform. Now, we give some of the N-Transforms and 
the conversion to Sumudu and Laplace [7,12]. 

Table 1: Special N-Transforms and the Conversion to Sumudu and Laplace 


fit) 

N|/(f)| 

s|/(r)l 

AM 

1 

1 

s 

l 

1 

s 

t 

u 

7 

u 

1 

7 

at 

e 

l 

1 

1 

s - au 

1 - au 

s - a 

, n— 1 

1 11 ~ 1 0 

n— 1 
U 

u"~ x 

1 

in-\Y n 

s n 

s n 

cos t 

s 

1 

s 

2 . 2 

S + U 

l + u 2 

1 + 7 

sinf 

u 

u 

1 

2 . 2 

S + U 

l + u 2 

1 + 7 


Some basic properties of the N-Transform are given as follows [7, 12]: 

If R(s ,u) is the natural transform and F(s)is the Laplace transform of the function f (t) , then 


N + [f{t)] = R(s,u) = - f f(t)dt = -f(- 
u ,, u \u J 


If R[s ,u) is the natural transform and G(u) is the Sumudu transform of the function fit), then 


N + [fit)] = Ris,u) = -]e 7 


f Ut ^ 


V S 


dt = -G 
s 


i„ 7 




If N + [/(()] = Ris,u), then N + [/(af)] = — Ris,u). 


If N + [/ (()] = Ris, u ) , then N + [/ it)] = -Ris,u)~ 


fio) 
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• If N + [/M]=i?M,then N + [/10] = 4M^)-4/(0)-^. 

u u “ u 

• Linearity property [12]: If Cl and b are non-zero constants, and / (f) and g{t ) are functions, then 

N + [af{t) ± bg(t)\ = aN + [f(t)\ ± b N + [g(t)] = a¥ + (s,u)± bG + {s,u). 

Moreover, F + (,s,w)andG + (,S,w) are the N-transforms of / (t ) and g (t ) respectively. 

4. THE NATURAL DECOMPOSITION METHOD 

The applicability of the natural decomposition method to Klein Gordon Equation is illustrated as follows: 

Consider the following Klein Gordon Equations (1.1) 

L t u(x,t) — L x u{x,t) + au(x,t) + F(u(x,t)) = h(x,t) 


Subject to initial conditions 

u(x,0) = f[x) ; w,(x,0)=g(x) 


(4.1) 


(4.2) 


02 02 

where L = — — ,L = — , Fu{x,t) is a nonlinear function of u(x,t) 3L!\d h{x,t) is a known analytic 

dt ' dx~ 


function. 


By taking the N-transform of Equation (4.1), we have 

N + [L f w(jc,f)]- N [L t w (.v', f )] + rr N ~ [u (.v, / )] + N 1 [Z 7 (u {x, t ))] = N [/; (.v, t )] 

Using the properties in Table 1 and the basic properties of the N-transforms, we get 
sw(x,0) 


— R(x,s,u)~ 


- N + [L v [w (x, t )]] + a N + [u (x, t )] + N + [/ (u t ))] = N + [h (x, t )] 


(4.3) 


substituting Equation (4.2) into Equation (4.3) to get 

sf{x) g{x) 


-R(x,s,u)- 


N + [L x [u (x, t )]] + a N + [m(x, t )] + N + [f(u(x, t ))] = N + [h(x, t )] 


sf{x) , g(x) 


—R(x,s,u) = — ^ + N + [L t [// (x, r )]] + N + [/; (x, t )] - a N + [u (x, r )] - N + [/ (u (x, t ))] 


R(x,s,u)= + + N + [L t [ M (x,/)]] + M 2 2 N + [/r(x,f)]- aM 2 2 N + [ M (x,/)]- N + [/( M (x,/))] 

s s s s s s 

Taking the inverse natural transform of Equation (4.6), we have 


(4.4) 


(4.5) 


(4.6) 
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N“‘ /?(a,5,m)=n" 


fix) ugix ) m 2 . T+ r, / 


2 2 

s s s 


+ N _1 


N + [L r [ M (x,0]] 


— N _1 


2 2 

W ^ T +r f w U 


a — N + [u (x, t )] + — N + [/ (u (x, t ))] 


(4.7) 


From Equation (4.7), we have 


ii(x,t)- j(x,t)+ N 


N + [L a (w(x,0)] 


■N 1 [/cm(x, /)] 


(4.8) 


u, z u. z 

Where ku{x,t)=Cl — — N + (u(x, I )) H — N [ / (u (xJ ))] and J (x.l) represents the term arising from the 


S S 

known analytical function and the given initial conditions. 

Now to deal with the nonlinear term, we represent the solution in an infinite series form. 


i{x,t)=Y J u n {x,t) 


n—0 


also, the nonlinear term 


ku{x, t ) can be written as ku{x,t)= ^ A n 


Where the A n 's are the polynomials of u 0 ,U { , ... , U n and can be calculated by the formula [4] 


A n = — — 

n\ dx" 


f\ 


n = 0,1,2, 


Substituting Equation (4.9) and Equation (4.10) into Equation (4.8) to get 


(4.9) 


(4.10) 


(4.11) 


^w„(x,f)= j(x,t)+ N 1 


-N H 




-2X 


(4.12) 


by comparing both sides of Equation (4.2) we conclude that 
U 0 {X,t) = j(x,t ) 


U x {x,t) = N 1 


u 2 (x,t) = N 1 


N + {L x u 0 (x,t)-A 0 ) 


l ^—N + (L x u ] {x,t))-A ] 


Continuing in this manner, we get the general recursive relation given by: 
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i{x,t) = N 1 


^rN + {L x u n {x,t))-A n 


, n > 1 


(4.13) 


Hence, from the general recursive relation in Equation (4.13), we can easily compute the remaining components 

u(x,t)as u 1 (x,t),u 2 (x,t),..., where u 0 (x,t) 


of 

given by 


is always the given initial conditions. Finally the exact solution is 


i{x,t)=Y J u n {x,t) 


n—0 


5. APPLICATIONS 


Here, we employ the NDM to three numerical examples and then compare our solutions to existing exact 


solutions. 


Example 5.1: Consider the homogenous Klein Gordon Equation 


«rr - - « = 0 


Subject to the initial conditions 

u(x,0) = 1 + sin x , u t (x, 0) = 0 

We first take the N-Transform of Equation (5.1), to obtain 

N ' [u n ] - N 1 [u xx ] - N 1 [t<] = N 1 [()] 

Using the properties in Table 1 and properties of the N-Transform, we have 

su{x, 0) u'(x, 0) 


-R(x,s,u)- 


N + kJ-N + [ M ] = 0 


Substituting Equation (5.2) into Equation (5.3), we have 


R(X, S,U ) : 


(l + sin x) 


+T N+ k+«] 


Now, taking the inverse N-Transform of Equation (5.4), we have 


(5.1) 


(5.2) 


(5.3) 


(5.4) 
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[/?(x,5,l<)]= N 1 

1 + sin x 

+ N _l 

T N + k*+“] 


s 


s 


(5.5) 


From Table 1, Equation (5.5) becomes 


<(A,r) = (l + sin;t)+N 1 


T N+ k + «] 


(5.6) 


From Equation (5.6), we can write 


<(x,f) = l + sinx+N 1 


— N + 
2 
S 


2X+IX 

n= 0 n = 0 


(5.7) 


Now from Equation (5.7), we can conclude that 


,(x,f) = 1 + sinx, M,(x,f) = N 1 


^N + [fl 0 +A 0 ] 


We continue in this manner to get the general recursive relation 


"„,iF')= n 


\n*[b„+a,] 


, n> 1 . 


(5.8) 


M, (x,t) = N 1 


"tN + [5 0 + A 0 ] 


= N 


-l 


^V N+ k« +« 0 ] 


Note that 


N" 1 

4 Nil] 

= N“‘ 

2 

U 

2 

rn 

= N _1 

2 

U 

3 


s 


s 

UJ 


s 


= — t 
2 


, (x,t) — 0 , V/t > 1. Hence u(x,t) = u Q (x,t)+ u l (x,t)+ u 2 (x,t)+ ... = 1 + sinx + — t 2 


and that 


u{x,t) = sin x + cosh t . 

Hence the exact solution is [1]. 

Example 5.2: Consider the inhomogenous nonlinear Klein Gordon Equation 


U n - U xx + u 1 = — X COS t + X^COSH 


(5.9) 


Subject to the initial conditions 


Index Copernicus Value: 3.0 - Articles can be sent to editor@impactjournals.us 



Natural Decomposition Method for Solving Klein Gordon Equations 


67 


u(x,0) = x , u t (x, 0) = 0 

We first take the N-Transform of Equation (5.1), to obtain 

N + [w„]~ N + [w lvr ] + N + [m 2 ]= N + [- xcost + x 2 cos 2 ?] 

Using the properties in Table 1 and properties of the N-Transform, we have 

su(x, 0) u\x, 0) 


—rR(x,s,u)-- 


u 


u 


u 


N + [» X J + N + [tt 2 ] = N + [- a:cos ( ] + N + [x 2 cos 2 ( ] 


Substituting Equation (5.10) into Equation (5. 1 1), we have 

x u 2 xf s ) u 2 x 2 ( s V u 2 + r 2 1 

+— N + k,-« 2 J 


R(x,s,u) = --‘ 


s s 


2 2 

ys +u~ j 


+- 


2 . 2 

\S +U 


Now, taking the inverse N-Transform of Equation (5.12), we have 


N 1 [/?(x, N’, r/)] = N -1 


- 


2 

U X 


s[s 2 + u 2 ) 

From Table 1, Equation (5.13) becomes 


+ N~‘ 


2 2 

U X 


_{s 2 + u 2 ) 2 


+ 


N + k-« 2 ] 


u{x,t) = xcost + x 1 sin 2 1 + N 1 

From Equation (5.14), we can write 

u(x,t) = xcost + x 2 sin 2 1 + N _1 




4n- 

s 


A, 

n = 0 n = 0 


(5.10) 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


(5.15) 

Here, A„ is the Adomian polynomial which represent the nonlinear terms. So, we compute few components of A„ 


and some values of B n . 


A q — u 0 

B 0 — u o„ 

A, = 2 u 0 u l 

5 1 = 

2 

A 2 = 2 u Q u 2 + m, 

^2 = U 2xx 


Now from Equation (5.15), we can conclude that 
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)(*,?) = xcost + A 2 sin 2 t , u x ( x , t ) = N~ 


^N + [fl 0 -A 0 ] 


We continue in this manner to get the general recursive relation 


«, ,,(.»,()= N 1 


An*[b„-aJ 


n > 1. 


(5.16) 


Note that we can calculate 


U { ( x , t ) = N 1 


^n*[s„-a] 


= N 


-i 


« XT+f 2 

— N Kx ,-«0 
s " 


= N 


-i 


2 

k“N + (2sin 2 1 - x 1 cos 2 ( -2a 3 cos t sin 2 1 + x 4 sin 4 t ) 

= N _I 

k“N + [- A 2 cos 2 1\ 

s 


s 


+ ... 


: ~ X ~ N“‘ 


f a 2 

u 


2 . 2 

\S +U 


+ ... = —x 1 sin 2 1 + ... 


u 0 ( x , t ) and u x { x , t ) 

Hence by cancelling the noise term that appears between one can see that the non-cancelled 

u 0 { x , t ) 

term of still satisfies the given differential equation, which lead to an exact solution of the form 

u ( x , t ) = XCOS t . 

This is in agreement with the result obtained by RDTM [2], 

Example 5.3: Consider the nonlinear, nonhomogenous Klein Gordon Equation 

2 t„ 2 o,2 . 4.4 


u tt — u xx + u — 2 x - 2 t +xt 


(5.17) 


Subject to the initial conditions 

u ( x , 0 ) = u t ( x , 0 ) = 0 


(5.18) 


We first take the N-Transform of Equation (5.17), to obtain 

N + k ] - N + [^ ] + N + [m 2 ] = N + [2a 2 ] - N + [2f 2 ] + N + [a 4 ( 4 ] 


Using the properties in Table 1 and properties of the N-Transform, we have 
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-R(x,s,u)- 


su(x, 0) u'(x, 0) 


' N+ [«„]+ N + [u 2 ]= N + [2x 2 ]- N + [2r ]+ N + [x V ] 


(5.19) 


Substituting Equation (5.18) into Equation (5.19), we have 

R(x,s,u)=2x 2 1 ^-j-4 1 ^y + 24x 4 + -u 2 ] 


s 3 s 5 


7 2 

5 5 


(5.20) 


Now taking the inverse N-Transform of Equation (5.20), we have 


[r(x,s,u)]=2x 2 N 1 

r 2 1 

U 

3 

-4 NT 1 

1 

£ 

_/i 4^. 

1 

+ 24jc 4 NT 1 

u 6 

7 

+ N _I 

t n+ L-« 2 ] 


5 


s 


s 


5 


(5.21) 


From Table 1, Equation (5.21) becomes 


u(x,t) = x 2 t 2 — — t 4 +— x 4 t 6 + N~ 
6 30 


"N + [m_„ -m 2 ] 


(5.22) 


From Equation (5.22), we can write 


u(x,t) = x 2 t 2 - — t 4 + — x 4 t 6 + N 
6 30 


-1 


— N h 

2 

S 


. n=0 n = 0 


(5.23) 


Here, A n is the Adomian polynomial which represent the nonlinear terms. So, we compute few components of A„ 
and some values of B„. 


A q — u Q 


Bq — u 0xx 


A, = 2 u Q u x B { = m 1x 

A 2 =2u 0 u 2 +u l 2 B 2 =u 2x 


Now from Equation (5.23), we can conclude that 


u Q (x,t)= x 2 t 9 - — t 4 +— x 4 t 6 , U l (x,t)= N 
6 30 


-i 


^n*[s„-a 0 ] 


We continue in this manner to get the general recursive relation 






n > 1. 


(5.24) 
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Note that we can calculate 


U l (x,t)= N 1 


4n1b 0 -a 0 ] 

= N _1 

« 2 m+ [ 2 ] 



5 " 


N" 


S 


0-2 . 2 2.6 4.4 . 2 2.6 1 6.8 . 1 4.10 1.8 1 8.12 

2 1 H — X t -x t H — x t X t H x t 1 X t 


= N 


-1 


^N + (2 t 2 -xY] 


+ ... = N 


-1 


15 90 

4^ -24a 4 4 

s s 


36 900 


, 1.4 1 4.6 . 

+ ... — — t X t + . 


6 30 


Hence by cancelling the noise term that appears between 

u 0 (x,t) and u^xj) u 0 (x,t) 

one can see that the non-cancelled term of still satisfies the given differential 

equation, which lead to an exact solution of the form 
u(x,t ) = X 2 t 2 . 

This is in agreement with the result obtained by MDM [1]. 

6. CONCLUSIONS 

In this research paper, the Natural Decomposition Method (NDM) was applied to solve three nonlinear Klein 
Gordon Equations. Exact solutions of the three applications were obtained. The method demonstrates significant 
improvement over existing techniques. 
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